Surface plasmon for graphene in the Dirac equation model 
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We consider single-layer plane graphene with electronic excitations described by the Dirac equa- 
tion. Using a known representation of the polarization tensor in terms of the spinor loop we show 
the existence of surface modes, i.e., of undamped in time excitations of the electromagnetic field, 
propagating along the graphene. These show up in the TE polarization and exist at zero tempera- 
ture. 



I. INTRODUCTION 

Surface plasmons are excitations of the electromag- 
netic field traveling along the interface between two me- 
dia or along a thin sheet, decreasing exponentially fast to 
both sides of the interface. These excitations have been 
studied in detail, both theoretically and experimentally, 
see, for instance [H, 0|- I n the present paper we show the 
existence of such excitations on graphene. The latter en- 
joys at present much attention because of its quite special 
properties, see, for example, The electronic excita- 
tions, primarily responsible for the interaction with the 
electromagnetic field, are best described by a Dirac equa- 
tion with correspondingly chosen parameters [|J. This 
holds, at least, near the Dirac points. The response of 
this model to electric and magnetic fields was convinc- 
ingly investigated in For a recent introduction to 
this model see @. 

Using a random-phase approximation with this model, 
surface modes on graphene were investigated in [7]. It 
was found that the existence of such modes is due to 
thermal excitations and that these do not exist at T = 0. 
In a number of papers this result was refined including, 
for example, electron-electron interaction [§| . In [t| [l(| 
surface plasmons were found in a conductivity approach. 

In general, surface modes exist on the interface be- 
tween two dielectric media with permittivities s± and £2, 
in case the relations £162 < and £1 + £2 < hold simul- 
taneously [ll| . For the surface of metal, described by the 
plasma model with 



(1) 



where lo p is the plasma frequency, these relations are sat- 
isfied for u < ujp/2. Such modes are known also on an 
infinitely thin metal sheet described by plasma confined 
to a plane, which may serve as a hydrodynamic model 
for graphene In this case, the response to the elec- 
tromagnetic field is described by matching conditions the 
field has to satisfy across the sheet. These conditions are 



quite simple and result in the reflection coefficients 
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(the superscript 'hy' denotes the hydrodynamic model, 
c is the speed of light), for a standard scattering setup. 
Thereby the electromagnetic field is separated into the 
usual polarizations, TE and TM, with plane wave ampli- 
tudes &(t,r) ~ exp(— iuot + ik||X||)$(s), 



$(z) 



e iqz +re -iqz } ( z < 0), 

te lqz , (z > 0). 



(3) 



Here the z-axis is perpendicular to the sheet, k|| is the 
in-plane wave number and q is the wave number perpen- 
dicular to the sheet. The Maxwell equations give 



lu 2 = c 2 (k 2 +q 2 ) 



(4) 



with fen = |k||| and the matching conditions result in the 
reflection coefficients ([2|) (and in corresponding expres- 
sions for the transmission coefficient t). For comparison 
we note the corresponding formulas for the case of dielec- 
tric media on both sides of the interface. In that case, 
the z-dependence of the amplitude is ~ exp(±igi j 2z) with 
wave numbers q\ and (72 in the respective media having 
permittivities £1 and £2. The Maxwell equations give 
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and the reflection coefficients arc 
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(5) 



A surface plasmon appears for real frequency ui, real wave 
number ki and imaginary wave number q, 



(6) 
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(or, correspondingly, imaginary 91,2), with real r\ which is 
the decay length of the amplitude in direction perpendic- 
ular to the surface, for a frequency giving the reflection 
coefficient a pole. 

In such setup, the spectrum of the electromagnetic field 
consists of scattering modes, having real all, lo, k|| and q, 
and the surface modes, in case these exist. For metals de- 
scribed by the plasma model and for graphene described 
by the hydrodynamic model, surface modes exist in the 
TM polarization. 
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It must be mentioned that the usage of the terminus 
'surface plasmon' is not unique. Frequently these are also 
called surface polaritions and excitations with complex 
frequency u), which are, strictly speaking, resonances, are 
included too. We mention also that surface plasmons 
need for their existence at least one translational invari- 
ant direction, i.e., at least the surface of a cylinder. If the 
surface does not have such direction, the poles of the re- 
flection coefficient appear at complex frequency and the 
mode decays. It must be mentioned that such modes can 
be quite long living and are of interest for applications. 

From the theoretical point of view, surface plasmons, 
as being part of the spectrum, are interesting per se. Spe- 
cial interest comes from the role of the surface plasmons 
in the Lifshitz formula describing the Casimir and van 
der Waals interaction, see [H, [lj| in the non retarded 
and [TH in the retarded cases for metals described by 
the plasma model and [l4j for the hydrodynamic model. 
Also these play an important role in the temperature cor- 
rections to the Lifshitz formula . Needless to mention 
their role in applications as nanophotonics. 

In the next section we show the existence of such plas- 
mons on graphene at zero temperature. 
Further in this paper we use units with h — c — 1. 



at the Minkowskian '3'-momentum 
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which is the relativistic invariant in the metric . These 
momenta are related by a Wick rotation, po = ip±, and p 
is the spatial momentum. In ([9]), a — e 2 /(47r) is the fine 
structure constant and N = 4 is the number of fermion 
species. The mass m is the gap parameter. 

The above momenta are related to the photon mo- 
menta in ([3]) and (0} by 



Po = w, 



E|| =p. 



Thus 



ui 2 — v 2 k 2 ^ 



(13) 



(14) 



will be real for tu > vkii and the relations 



P\\ = \ + K =r ) 
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hold with real r\ and p implying the restriction uj < fen 
Using these notations, we rewrite 



II. SURFACE PLASMONS IN THE DIRAC 
MODEL 

The description of graphene by a Dirac equation is a 
model for the electronic excitations accounting for the 
linearity of their spectrum near the Fermi points. In the 
language of quantum field theory one has to calculate a 
spinor loop in (2+1) dimensions in a metric 



g^ v =diag(l,-i>,-u), 



(7) 



where v = 1/300 is the Fermi speed (in units of the speed 
of light), and to couple the emerging polarization tensor 
to the electromagnetic field in (3+1) dimensions. This 
was done in a number of papers. We use the representa- 
tion in [Hi at T = (it was extended in [18] to T ^ 0), 
since here the expressions for the reflection coefficients in 
terms of the polarization tensor were derived. At T = 0, 
these reflection coefficients arc 
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Qtm = P^3(p), 



(8) 
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with the Euclidean momenta 



?TE = ^*(P). 



At this point the signs become important. First, we note 

$(p) < (17) 

for < p < 2m, which can be seen explicitly from (|lip . 
It should be mentioned that this condition, together with 
the reality of 77, defines the only region where Qte and 
Qtm are real. The upper bound is the threshold of pair 
production. We emphasize that &(p) takes negative val- 
ues. This is due to the extra minus sign a spinor loop has 
from the statistics, in distinction from a bosonic loop. 
This sign makes the reflection coefficient of the TE po- 
larization having a pole, and as a consequence, a surface 
plasmon appears in the spectrum. 

For comparison we mention the hydrodynamic model. 
Using the same notations as in (fit?)) we have 
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and the pole is possible only in the TM polarization. 

Thus, in the Dirac model, we may have a pole in the 
reflection coefficient of the TE polarization. Its location 
is solution of the equation 
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Inserting from (|16p and (jlll) , this equation can be rewrit- 
ten in the form 
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FIG. 1: The solutions of eq. ((20)) for v = 1/2 and several 
values of the parameter a. From bottom to top, the curves 
correspond to a = 1, 0.5, 0.2, 0.1, 1/137. 
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FIG. 2: The solutions of eq. (J20j) for v = 1/300 and several 
values of the parameter a. From bottom to top, the curves 
correspond to a = 1, 0.5, 0.2, 0.1, 1/137. 



where we used the first line in eq. (|15l) for r\ and have 
from the second line 
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v 2 k 2 
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This equation defines a relation between ui and fen. Its 
solution, 



(22) 



is the dispersion relation of the surface plasmon. For this 
solution, from the reality of r\ and p and using eq. (|15|) , 
the restriction 



vkii < oj s f(fc||) < fen 



(23) 



follows. Further, the solutions does not exceed the 
threshold of pair creation, i.e., they satisfy 



p < 2m, 



(24) 



since otherwise $(p) would not be real and the equation 
([20)) could not have a real solution. 

The solution of eq. ([20)) exists for all fen. It can be 
found numerically. Such solutions are shown in Fig. 1 
for v — 1/2 and in Fig 2 for v = 1/300 for several values 
of a. 

The asymptotic behavior for small fc| | can be found by 
rewriting eq. (|20[) in the form 
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and iterating starting from inserting u) = fcn in the right 
hand side. One obtains 



w«f(fc||) - fe|| - r )- - fc J 
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which is the expansion of the solution in powers of k\ \ . 
In order to get the behavior for large feii, it is meaningful 
to rewrite the equation in the form 
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v 2 kf, + (2m) 2 
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(27) 



It can be iterated by inserting ui = ^Jv 2 k 2 ^ + (2m) 2 in 
the right hand side. The solution is 

Wsf (fc||) = ^ 2 fc 2 + (2m) 2 + ... , (28) 

where the dots denote contributions exponentially small 
for large fen. 

In this way, we have for both limiting cases a linear 



dispersion relation, 

{fcii for ->■ 0, 
(29) 
v.k\\ for — > oo. 

The behavior in between these limiting cases, as can be 
seen from the figures, is smooth for any fixed values of the 
parameters a and v. For small a, the curve has a quite 
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sharp knee. For the physical values of the parameters 
it can be considered as consisting of two straight lines, 
crossing in Am = 2m. 

The massless case corresponds to large both, uj and k\ \ . 
In that case we get from eq. ([28]) 

w s f(fe||)| m= o = vk\\, (30) 

i.e., a linear dispersion relation. It must be mentioned 
that this limit cannot be performed in the function &(p), 
resulting from the polarization tensor for any finite p, 
because, as we know by hindsight, p becomes small not 
to exceed the pair creation threshold ([24l . 

III. CONCLUSIONS 

We have shown the existence of a surface plasmon on 
graphene. We used the Dirac model, accounting for the 
spinor loop without further approximation. The surface 
plasmon shows up in the TE polarization which is due to 
the minus sign a Fermi loop has as compared to a bosonic 
one. 

The surface plasmon appears at a frequency giving the 
reflection coefficient a pole. The corresponding equation, 
(|20[) . defines the frequency cj s f(/c||) as a unique function 
of the in-plane wave number fcii. This is the dispersion 
relation for this plasmon. For small and large fen, eq. 
(|29|) . and in the massless case, eq. (|30|) . this relation is 
linear. 



We found the surface plasmon at zero temperature, 
T = 0. It is clear that it survives also at finite tempera- 
ture if that is sufficiently small. This is, because tempe- 
rature corrections to Qte (as well as to Qtm), are always 
positive. These are small at sufficiently low temperature, 
but grow with increasing temperature until exceeding the 
negative Qte- 

We have to mention that the plasmons found here are 
different from those known in literature. The closest to 
ours are those in [10] which exist in the TE polarization. 
However, the frequency condition used there (eq. (1) in 
[nil ]) looks somehow oversimplified as compared with eq. 

Finally we mention that the surface plasmon found 
here, has a good chance to exist on a carbon nano tube 
too. In the present paper we assumed a flat graphene 
sheet, neglecting the ripples. We did not estimate their 
influence. Also we did not discuss any interaction be- 
tween the electrons in the graphene, or any interactions 
with the lattice which are beyond the considered model. 
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